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Abstract 



We measured the two dimensional galaxy-galaxy correlation function from the Muenster 
Red Sky Survey (Ungruhe 1999). This survey has a slightly larger surface as the APM 
survey and complete up to rp <18.3 mag (z me d ~0.14). The large dynamical range of the 
survey made it possible to examine sub-catalogues with different limiting magnitudes. The 
positive part of our measured angular correlation functions (after the colour correction) 
agree excellently with those measured by Maddox, Efstathiou and Shutherland 1996 from 
the APM galaxy catalogue (Maddox et al. 1990). From the measured angular correlation 
functions the three dimensional two-point correlation function £(r) was calculated. We 
have chosen the evolution parameter e so that the change of the correlation length ro with 
the depth of the sub-catalogue were minimal. This prescription gives ro = 5.82±0.05 h _1 
Mpc and a very high, positive evolution parameter e = 2.6. The position of the point, 
where the 3d correlation function breaks away from the power law depends on the depth 
of the sub-catalogues: for low magnitude limits the break-down comes earlier, for higher 
later. 



1. Introduction 

Since almost half a century the two-point correlation function has been the most popular 
statistic to describe the galaxy clustering in astrophysics. It has been firmly established 
(Peebles 1980 and references therein) that the three dimensional two-point correlation 
function (3d CF) is of the form £(r) = (ro/r) 7 in a wide range of distances, breaks away 
from the power law form at certain distance r\ and becomes negative at 7*2. The 3d CF can 
be measured by using magnitude or volume limited three dimensional galaxy catalogues or 
can be calculated from the 2d CF, measured from large surface two dimensional catalogues. 
The parameters ro and 7 measured from local and medium deep catalogues are ro ~ 5/i _1 
Mpc 7 = 1.7 — 1.8. Some authors estimate also the parameter T2 ■ Baugh 1996 finds by 
using the APM galaxy survey (Maddox et al. 1990) r 2 « 30/i _1 Mpc, Tucker et al. 1997 
measure from the LCRS (Shechtman et al. 1996) T2 ~ 30 — 40/i -1 Mpc while Cappi et al. 
1998 on the basis of the Southern Sky Redshift Survey (SSRS2, da Costa et al. 1994) find 
that the first zero point of the autocorrelation function £(r) increases linearly with the 
sample depth. Several authors notices that the 3d CF rises above the power law before it 
breaks down (Baugh 1996, Guzzo et al. 1991). 

To measure the detailed form of the tq{z) function large surface deep 3d catalogues are 
needed. At present the deep 3d surveys (say z me ^ > 0.3) contain several hundred galaxies 
and are not large enough for this purposes. Usually the evolution of the CF is parametrised 
and the parameter is measured. The customary parametrisation of the correlation length 
ro as a function of the redshift is 

r (z) = r Q (l + z) = T*. (1) 

This relation was proposed in 1977 (Groth & Peebles 1977) when the observable redshifts 
were much smaller than 1. For z <C lthe first term of expansion of tq(z) by z is Tq(z) = 
ro(0)(l — ^-^z) (ro measured in physical coordinates). However, for large z with e = const. 
(|l]) has no solid physical basis. In a model proposed by Ma 1999 e varies from 7 — 1 at high 
redshift trough at medium z down to 7 - 3 for low redshift. The first value corresponds to 
the linear grows of the perturbations (over the general expansion) , e = if the structures 
are fix in comoving coordinates and e = 7 - 3 when the clustering is stabile in physical 
coordinates. The canonical value of the correlation length is ro(0) = 5.4 h Mpc (Davis & 
Peebles 1983). From small catalogues in general, one can only measure only a combination 
of ro and e (Le Fevre et al. 1996, Carlberg et al. 1997, Connolly et al 1998, Small et al. 
1999) therefore it is difficult to compare the results. The e evolution parameter can be 

roughly estimated from the equation ro(z m ed) = r o(0)(l + z me d) 1 , where z me( i is the 
median redshift of the catalogue, ro measured in physical coordinates. It is clear that the 
amplitude of the 3d CF decreases strongly with increasing redshift and for fix ro(0) the 
measured e is between and 2.6. Postman et al. 1998 determined the 3d CF from an I 
selected 2d catalogue. From the I<20 sub-catalogue they found ro = 5.2 ± 0.4/i _1 Mpc 
and e positive. Deeper the catalogue (I<23, z me d = 0.5), e is smaller ie. the evolution 
is slower what is in contradiction with the results of Ma 1999. Small et al. 1999 find no 
evidence that the comoving correlation length changes between the redshift intervals (0.2 
- 0.3) and (0.3 - 0.5). 

From a large surface catalogue with large dynamical range (as APM or MRSS) one can 
separate a number of sub-catalogues eg. different limiting magnitudes m. In this case by 
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taking the same ro for every sub-catalogue, one can determine the parameter e (at least 
for z <Cl). 

A further important question is, where the CF breaks away from scaling, at which 
distance r becomes negative and are there measurable anti-correlations and correlations 
after that point. These problems are related with the question of the size of the largest 
structures in the Universe. Wu, Lahav & Rees 1998 and 1999 find using several indicators 
that on scales larger than 300 h _1 Mpc the fluctuations become negligible. 

Pietronero and coworkers (Coleman & Pietronero 1992, Sylos-Labini et al. 1998, Joyce 
et al. 2000 and references therein) question the applicability of the CF to characterise the 
large scale structure of the Universe. They state, on the basis of conditional probability 
analysis of different catalogues that the correlation length ro for clustering of galaxies do 
not characterises intrinsic properties of the distribution, but the particular sample one is 
considering and no statistical evidence for the usually assumed homogeneity on a scale 
large enough. 

In this paper we determine the two-point correlation function from the Muenster Red 
Sky Survey (MRSS, Ungruhe 1999). Section 2 contains a brief description of the survey. 
Section 3 presents the measuring methods of the 2d CF, the analysis of the results and 
comparison with the APM 2d correlation function as determined by Maddox, Efstathiou 
and Shutherland 1996 ( in the following MES). In section 4 we calculate the 3d CF from 
different sub-catalogues using the Limber-equation, discuss the scaling properties of the 
CFs and determine the correlation length ro, as well as the evolution parameter e. The 
interpretation of the results is discussed in Section 5. 

The two-point CFs give only partial information about the galaxy distribution, for 
more detailed description higher order CFs are needed. On this correlation function we 
report elsewhere. 

2. The catalogue 

The Muenster red sky survey (MRSS, R. Ungruhe 1999) is based on 217 contiguous 
plate of the ESO Southern Sky Atlas R, covering more than 5000 square degree. The 
plates were digitalised using two PDS 2020 GM phs machine of the Astronomical Institute 
Muenster. The image surface brightness profiles were used to distinguish galaxies from 
stars and other objects. The catalogue contains more than 7 million galaxies and out 
of these about 2.3 million brighter than rp=19.5 mag, but to avoid the "last magnitude 
effect" we do not use in the determination of the CF the galaxies fainter then 18.5 mag. 
After corrections the magnitudes are accurate to < 0.05 mag across each of the 217 plate. 
Images in the overlapping regions of neighbouring plates are used to establish a uniform 
magnitude system over the entire survey. In order to calibrate the matched magnitudes, 
CCD photometrie of 1766 galaxies and 1876 stars, homogeneously distributed on 92 plates 
were used. The estimated plate to plate zero-point-error is less than 0.10 mag. The 
extinction was corrected using the maps of Schlegel et al. 1998. The catalogue is complete 
up to 18.3 mag (Ungruhe 1999). 

The catalogue is complemented with galaxy images digitalised in pixels with 1.01 arcsec 
edge length. The pictures of objects were stored in frames (out-cuts) with side length 21, 
51 and 101 pixels, choosed by the software. The computer classification of about 2.6 
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million objects were carefully controlled visually (Ungruhe 1999), what reduce essentially 
the possible star-galaxy misclassification. 

The corresponding J-atlas was digitalised at two UK institutions: at the Institute for 
Astronomy in Cambridge using the APM machine (186 Schmidt-plate, Maddox et al. 
1990) and at the Royal Observatory Edinburgh (154 Schmidt-plate) with the COSMOS 
machine Stevenson et al. 1985. Both catalogues are included in MRSS that makes possible 
the comparison of the positions and other quantities, as well as the determination of the 
rp - bj colours of the galaxies. Because of the large surface of the survey the integral 
corrections are negligible. 

The F band has some advantages over the J band: the vignetting and desensitization 
corrections are 2-3 times smaller, consequently the magnitudes towards the edges of 
plates are better determined. Furthermore, the K - correction is also smaller and varies 
mildly with the depth of the catalogue. In spectral range of the MRSS the average rp - bj 
is 1.5 mag (Jones et al. 1991). In addition the MRSS catalogue contains other data such 
as surface brightness, ellipticity, length of small and large axes, large-scale alignments etc. 
Schuecker et al. (1990, 1994) measured the objective prism redshifts for 900 000 galax- 
ies, and Spiekermann 1994 gave an automated morphology classification of all galaxies 
brighter than rp = 18.5 mag. 

3. Measuring the correlation function 

As it was emphasised by Maddox, Efstathiou and Shutherland (1996), large areas of 
the sky are required to determine the 2d CF w reliably on angular scales of a few degrees. 
They used 80 to 120 contiguous plates to compute the CF. Here we determine the w(9; m) 
from 152 contiguous plates. 

For measuring the galaxy - galaxy correlation function we use two essentially different, 
independent methods: (i) count in annuli around galaxies and (ii) count in cells. 

3.1 Count in annuli around galaxies 

The most direct method to measure the over-density around galaxies is to count the 
number of galaxies between circles of radius 0j and The corresponding density can 

be calculated by dividing by the surface of the annulus. If the first galaxy is near to the 
border of the sample only a part of the annulus is taken into account being inside the 
sample. For very small angular distances this method is more precise than measuring 
the surface by the number of random points in the annulus. This method does not need 
gridding of the data, it is simple and fast to apply but, as it was shown (Hamilton 1993, 
Maddox et al. 1996), subject of first order errors in the galaxy density contrast. We use 
this estimation only for small angular scales 6 < 1°, when the amplitude is high. 

3.2 Count in annuli around bins 

For larger angular distances we do not need the exact position of the galaxies. We 
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binned the galaxies in square cells in an equal area projection, and applied the count in 
annulus method for the cells. A cell belongs to a given ring when its centre is inside the 
ring. The surface of an annulus is taken to be equal to the the sum of the surfaces of cells 
in the annulus. That is we use the usual estimator 

= < n i n j > L (2) 

< Hi >< n,j > 

In the angular range 0.08° 9 < 0.8° we use a grid of edge length s = 0.01° , while for 
0.5° < 9 < 4° we have s = 0.1° and s = 0.25° for 9 > 4°. The overlapping parts of the 
CFs calculated in different intervals agree well with each others. 



3.3 The count-in-cells method 



Another possibility to determine the CF is to measure its average on square cells of 
edge 9: 

w ( ) = ^ / d2e ^ J d 2 02M|0i - e 2 \). (3) 

This quantity can be expressed in terms of the moments of galaxy-number in a cell as 
follows (Rubin 1957, Peebles 1980) 

-/m N*-N 2 -N 

w{9) = — 2 . (4) 



To compute iV = N(6) we put a grid with edge length of 9 on the equal-area-projection. 

In order to be sure not to cut structures in parts we shift the grid in steps <C 9 in 5 
and a directions and count the galaxies in cells for every position. By calculating the 
averages the catalogue is heavily oversampled. The - mostly beneficial - consequences of 
oversampling was discussed by Szapudi et al. 1997. 

From the average of the CF one can determine the true CF by inverting the integral (|3|). 
If the CF is of the form w (9) = A9 1 ' 1 then its average can be calculated by numerical 
integration from (|j|): 

w Q {9) = 1.968A0 1 - 7 = 1.968u/ o (0). (5) 



3.5. Variation of w{9) with the limiting magnitude 



By computing the CF from sub-catalogues with different limiting magnitudes we can 
mimic the changes in depth. The limiting magnitude rawas varied from m = rp(max) = 
15.5 mag up to m=18.5 mag in half magnitude steps. The result is shown in Figure 1. 
The CFs on a broad, on a log-log plot more than two decades range are straight lines, i.e. 
w(9,m) = A(m)9 1 ~" ( . The parameter in the exponent is 7 = 1.69 ±0.02. Using this value 
of 7, we can measure the amplitudes of the straight part of the CF. The parameters 
log (A(m)) deg 7 " 1 ) are given in the Table 1. 
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Table 1. 



m 


15.5 


16.0 


16.5 


17.0 


17.5 


18.0 


18.5 


\og{A deg 7 " 1 ) 


-.64 


-0.79 


-0.95 


-1.12 


-1.28 


-1.465 


-1.63 



In the 15.5 < m < 18.5 interval the amplitude of the CFs A(m), measured from the 
horizontal part of the function w(9;m)/9 1 ~^ can be parametrised as 

A(m) = 0.0148 • io - 33 ^ 15 - 5 -™^ 7 - 1 . (6) 

The amplitude changes a factor of ten in the three magnitude range, as in MES or Postman 
et al. 1998. The whole 2d CF, can be parametrised as 

w p (6;m) = A p (m)9 1 -^ 

We use a = 9 and v = 2.5. The first bracket on the right hand side describes the break- 
down of the scaling at small angles, while the break away from the power law at large 
angular distances contained in the second. In this formula we take into account that for 
small limiting magnitudes the 2d CFs are somewhat steeper, and 7 P varies slightly with 
the limiting magnitude, consequently A p (m) does not agree exactly with A{m) given 
above, where we have used the value 7 = 1.69. The parameters in (0) are given in Table 
2. 



1 + 



1 + 



1 + 



(7) 



Table 2. 



m 


A p deg> _1 


7 P 


6>i deg 


2 deg 


(9 3 deg 


15.5 


0.2352 


1.73 


8.8 10~ 3 


6.5 


40 


16.0 


0.1740 


1.71 


7.1 10~ 3 


8.5 


40 


16.5 


0.1350 


1.69 


5.5 10~ 3 


8.0 


40 


17.0 


0.0860 


1.68 


4.5 10~ 3 


7.5 


40 


17.5 


0.0550 


1.68 


4.1 10~ 3 


7.0 


40 


18.0 


0.0350 


1.68 


3.7 10~ 3 


6.3 


40 


18.5 


0.0245 


1.68 


3.4 10~ 3 


6.0 


40 



The maxima of the CFs are shifted toward smaller angles as the limiting magnitude 
grows. This corresponds to the change of the average apparent galaxy diameter with m. 
When the apparent picture of the galaxies overlap, the software identifies them in the 
"unknown objects". One can parametrise the maxima of the CFs as function of m as 
6» max = 2.83 • HTelk. 

The slope of the straight-line part of the log(u>)-log(#) function is between -1.73 and 
-1.68 in agreement with earlier results. The functions with lower limiting magnitudes 
are somewhat steeper. This could mean that the intrinsically brighter galaxies are more 
correlated (Hamilton 1988, Davis et al. 1988, Benoist et al. 1996). 

Between 6° and 10° the CFs go through zero (Figure 12 - 14). The amplitude of the 
negative part decreases with increasing magnitude limit. 
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3.6. Consistency of the two methods 

The CF measured by count in annuli method and the average of the CF measured by 
counts-in-cells for galaxies brighter than rp = 18 mag are shown in Figure 2. First of all 
we check whether the two methods are consistent by computing the average of the CF 
determined by count in annuli method. 

The broken line in Figure 2. shows the the parametrised CF w p (9). We compute the 
average of this function on squares with edge length 9. The direct Monte Carlo integration 
of w would be dangerous because of the peek at 9 = 0. The average of wq{9) is known 
and the average of the difference w p {9) — wq{9) is not singular any more and can be easily 
calculated from (|7]). The resulting averaged CF is shown by the dashed line in Fig. 2. 
The agreement of the two methods is excellent. 

3.7 Comparison with w(9,m) from the APM survey 

Maddox, Efstathiou and Sutherland 1996 have compared in detail the CFs, calculated 
from the APM survey with w(9, m) from other surveys. An important result of their 
analysis that large area surveys are required to determine w reliably on angular scales of a 
few degrees. They found that the deviation of w from a power law seen in the small area 
catalogues are caused by the integral constraint in the estimator of w. 

The APM CF agrees well with that calculated from the Lick-survey (Groth and Peebles 
1977) for 9 < 3°, but the APM CF breaks less sharply from the power law on scales larger 
than ~ 3°. By comparing with the CF calculated from the Edinburgh-Durham Southern 
Galaxy Catalogue (Collins et al. 1989, 1993) MES found that up to 2° the two CFs agree 
well but for larger angular distances the EDSGC gives more correlation than the APM 
catalogue. 

In Figure 3 we compare our CFs with that of MES (their Figure 23). The magnitude 
limit (after the colour correction bj — rp = 1.5 mag) for the two upper curves on both 
panel are not exactly the same and this causes the small relative shift of the CFs. The 
agreement of the other pairs of curves are excellent. We confirm the results of MES con- 
cerning the 9 and m dependence of the 2d CFs. Note that the APM CFs were calculated 
by using 80 - 120 Schmidt plates, the MRSS 152 plates; the plate sizes and their overlaps 
in the two catalogues are not the same, furthermore the N(bj) and N(rp) functions are 
different, the photometry and the star-galaxy separation software are independent etc. In 
spite of all these differences the CFs measured from the two catalogues agree in a range of 
several magnitudes. After the colour correction our CFs, determined from sub-catalogues 
with different limiting magnitude m, agrees excellently with the CFs measured from the 
APM catalogue. For lower magnitude limits (m < 17.0 mag) the MRSS functions show 
less correlation in the region 9 > 4°. MES give the negative part of the 2d CF only for 
one magnitude limit. Their CF shows much less anti-correlation as ours (Figure 4). 
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3.8. Scaling of the correlation functions 

Groth and Peebles 1977 have shown that the CFs of different mu m can be transformed 
into one curve by shifting them with appropriate amounts A x (m) and A y (m) in the log(0) 
and log(w) directions. As it is shown in Figure 5 one can approximate the measured CFs 
on log-log scale with straight lines with slopes 5\ = -0.68 and 62 = -3.0. We could scale 
the section point into one points, in accord with to the Groth - Peebles 1978 prescription. 
When scaling the m = 18 mag CF to m = 17 mag we must use A x = 0.0055, A y = 0.332 
and for m = 16 mag CF to m = 17 mag the parameters A x = —0.0038, A y = —0.231. 
For the APM catalog MES obtained a reasonable overlap of the transformed CFs. Using 
the EDSGC catalogue CNL92 demonstrated similar scaling. However, it will be shown by 
assuming the scaling of the CFs that the shifts in the log(#) and log(w) directions are not 
independent from each other. We discuss further the scaling properties in the section 4.2. 

An other way to demonstrate the scaling is to divide the measured CFs by A p (m)9 l ~' lv ^ . 
The resulting function f(9), is independent of the limiting magnitude m in a large angular 
range, at least up to 6° (Figure 6). The semilogarithmic presentation of the f{6) function 
emphasises the dip at about 0.4°. A similar dip is seen in Figure 23 of MES96 and Figure 
2 of CNL92. The reason for this dip seems to be obscure. 

The CFs measured from MRSS and from the APM catalogue deviate from the power 
law at about 2°, independent of the limiting magnitude. The CFs determined from EDSGC 
break away from scaling at about 19°, also independent of m (CNL92 Figure 2). 

The fact that the function f{6) does not depend on m in a bright range of 9 shows 
already that the 3d CF, calculated from the 2d CFs, must depend on the limiting magni- 
tude m, i.e. on the depth of the catalogue. 

4. Relation between two dimensional and three dimensional correlation functions 

From a large surface deep 3d catalogue one could determine the luminosity function 
of the galaxies, together with their luminosity evolution, as well as the 3d correlation 
function and its evolution with the redshift. The 2d catalogues, like MRSS, contain at 
least one order of magnitude more galaxies but the problems of the selection function is not 
solved. In some sense the informations got from the two type of catalogues complement 
each other. 

Here we determine from the measured 2d CF the ro and e parameters of the 3d CF. 
In the literature there are several way to simplify the problems. One possibility is to take 
a dN/dz function, measured from a smaller catalogue, fix the parameter e and calculate 
ro as a function of the magnitude limit. The second way is to fix ro and e and choose 
appropriate parameters in the selection function. A further method is to use a reliable 
dN/dz function and to choose e so that ro doesn't depend on the limiting magnitude 
(depth) of the sub-catalogues. 



7 



4.1 The Limber-equation 



The Limber-equation connects the 3d CF with the 2d CFs: 



e 

where = 2sin(0/2), and the kernel is 



w{@;m) = — I rK(r/Q;m)£,(r)dr (8) 



o 



_J_ /* F(x)(rfiV/^ & 

(iV^io (1 + *)*+«■ lj 

Here u> is the solid angle covered by the survey, N*(m) is the number of galaxies per 
steradian, F{x) — 1 is a relativistic correction (Peebles 1980, §56), which is exactly zero for 
Qq = 1, and small for any CIq for the depth covered in MRSS. For the simplicity we use 
the Einstein- de Sitter model (f2o = 1, A =0). In this case the metric distance depends on 
the redshift as 

x = ^i-VTT~z). (io) 

We use H = WOh km/s Mpc. 

If the 3d CF is of the form £(r) = (^) 7 then from (g) and © follows that w(9; m) = 
A(m)9 1 ~~ t . The amplitude A(m) can be calculated from £(r) and the normalised galaxy 
number distribution p(x;m): 

dz 1 dN(z;m) dz 
p(x;m) = q{z;m)— = —— — ; 11) 

dx N*(m)uj dz dx 

with 

roo coo 



p(x;m)dx = / q(z;m)dz = 1. (12) 
Jo 

The amplitude A(m) dependent only on the form of the ^ function and independent of 
its amplitude. If we know p(x;m), there is a unique connection between the correlation 
length ro and the 2d amplitudes A(m). 

4.2. Scaling and the Limber-equation 

Let us now assume that the 2d CFs scale with the limiting magnitude m, ie. 

w{0;m!) = Cw{\e-m). (13) 
In this case from (0) follows that K (t; m) scales as 

K(t-m') = jK(^;m) (14) 



and by ([ll]) we have 

A 



p(x; m') = \l —p(x/\; m) 



1 + z{x/X) 
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From the normalisation condition of p(x; m!) follows that 

Alc-5j(A,m,e) = 1, 

where 

/■oo m + z(xA)l (3+e)/2 
I(\,m,e)= — — — — p(x;m)dx. 
Jo 1 + zlx) _ 



(16) 



(17) 



That means that the shift of the w(6;m) functions in directions of the log(#) axis and 
log(u> ) axis are not independent but 



AlogH = log(C) = 3Alog(0) + 21og(/(A, m, e)). 



(18) 



For a given realistic distribution function p(x; m) or q(z; m) the integral /(A, m, e) can be 
calculated. It turns out that it gives only a small contribution to C, corresponding to the 
fact that in the depth of MRSS the curvature effects are neglected and the redshift effects 
are small (see Groth & Peebles 1977). 

Consequently, if we scale together on the log-log plot the less steep part of the CFs 
taking into account (|i~8|), their steeper parts will not fall together. That means that the 
correlation amplitudes and the position of the break in w, as measured from the MRSS 
are independent, ie log[w(6; m)\ can not be transformed into log [w{6\ m')] by a shift of 
Alog(6>) and Alog^) given by (^). We conclude that if there exists a 3d CF, independent 
of the depth of the catalogue, the 2d correlation functions do not scale with the limiting 
magnitude m. The slope of the 2d CFs and the position of the break down of the scaling 
are independent from each other. 



4.3 The selection function 



Maddox et al. 1990 proposed a selection function calculated from a redshift dependent 
luminosity function. Later, MES96 used a selection function deduced from the measured 
redshifts of a part of the galaxies in the 2d catalogue in the form of 

(IN 

— = N*ujq(z;a,(3,z ) (19) 
dz 

where q(z; a, (3, zq) was parametrised as 

q{z\a,P,zo)= , P . fi-Ve-'t' 8 , (20) 



r (^) 



20 



with a — 2 and (3 — 3/2. The maximum of this distribution is z max — (a/P^^zo, the 
mean value of the n-th moment of the redshift z n {m) is 

(2D 

r((a + l)//3) 

and the modal value can be calculated from the following equation 

r(y,u4 r (y>' (22) 
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For a, — 2 and (3 — 1.5 one has z max — 1.2114zo, z — 1.5048 zo and z moc j — 1.412;zo. From 
these dN/dz function one can determine the luminosity function. It depends mildly on z. 
This parametrisation is simple and general enough to describe the selection function. 

We use two set of parameters a, (3 and zq. First we take, following MES96 a = 2 (3 = 1.5 
and determine zq from the condition rg = 5/j -1 Mpc (and e = -1.3). The second set of 
parameters was calculated by using the luminosity function, determined from the Las 
Campanas survey (Lin et al. 1996). This must better approximate the real selection 
function, because the R band used in the Las Campanas survey is not far from the rp 
band of the MRSS. The parameters are given in Table 3. 

Table 3. Parameters of the selection function 





Parameter set 1 
tq = 5/i -1 Mpc 


APM 
bj < m+1.5 


Parameter set 2 
From Las Campanas LF 


m 


a 


P 


zo 


zo 


a 


P 


zo 


z 


15.5 


2 


1.5 


0.0266 


0.0326 


2 


1.50 


0.0284 


0.0436 


16.0 


2 


1.5 


0.0327 


0.0366 


2 


1.50 


0.0350 


0.0527 


16.5 


2 


1.5 


0.0414 


0.0439 


2 


1.52 


0.0440 


0.0651 


17.0 


2 


1.5 


0.0516 


0.0534 


2 


1.57 


0.0568 


0.0804 


17.5 


2 


1.5 


0.0645 


0.0646 


1.95 


1.62 


0.0730 


0.0965 


18.0 


2 


1.5 


0.0834 


0.0774 


1.90 


1.66 


0.0921 


0.1205 


18.5 


2 


1.5 


0.1050 


0.0915 


1.90 


1.695 


0.1140 


0.1461 



In the second parameter set the parameter (3 grows with the depth of the catalogue. This 
fits in the general picture. For deep {z > 0.3) catalogues several author (Lilly et al. 1995, 
Brainerd & Smail 1998) use the same function form with j3 = 2, or modifications of it( 
Postman et al. 1998). We note already here that our main results do not depend essen- 
tially on the exact form of the selection function. 

4.4. Solution of the Limber-equation for CF of power law form 



When the 3d CF follows a power law and the normalised selection function is known, 
ie. the parameters a, (3 and zq are given, there is a simple relation between A{m) and ro 
(Totsuji et Kihara 1969). For small r the 3d CF is £(i8) = (j§) 7 , and we get 



POO 

A(m) = @^- l w{Q- m)=rl / t^Kfa m)dt. 

Jo 



(23) 



By using ( |i~0| ) and developing ^§(l + z) ^ 3+<E - ) by z up to the second order, we can compute 
the integral (^3|) : 



A(m) = C 



where 



cz 
2 + 2a - 7 



r (r) + dr (r - 1/(3) z Q + C 2 r (r - 2/(3) z 2 Q 



rW)rm 



(24) 
(25) 
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and 



_ 1-9 + 37 
2 /3 ' - 



4e 



_l9 7 2 
2 f^— 



657 + 16e 2 + 88e - 24 7 e + 116 



32 



(26) 



With the selection function, calculated from the Las Campanas luminosity function (Lin 
et al. 1996) we can determine from the measured A(m) using (|24) the parameter rg: 



c 1 
r = — z A(m)-/C~ 
-no 



r(r) + dr(r - l//3)ab + C 2 r(r - 2//3)zg 



(27) 



If we fix the correlation length ro and e, the equation (^J) can be solved for zo by 
iterations. First we neglect the terms containing zq in the square bracket. Then 



(0) . , .-lH ro 
Zq = A(m) 1 [CT(r)J- 



(28) 



substituting that into (p3) one gets 



(m) 



l + Ci- 



,(0) 



T(r) 



r(r-2//3) (0) 2 

r(r) z ° 



(29) 



For ro = 5h^ 1 Mpc and e = -1.3 Zo(m) is shown in Table 3. and can be parametrised as 
zo(m) = 0.024 + 0.0117(m — 15. 5) 1 ' 63 . Of course, it is not granted that these zq values 
correspond to the real dN/dz distribution. This parametrisation agrees well with that of 
MES for m ~ 17 mag, but our zo(m) changes somewhat faster with m. The same is true 
for out second parametrisation. 

Next, with the second parameter set we compute ro(m), for e = —1.3. In Figure 7. 
we show our ro(m) function with and compare it with the same function, calculated from 
the data given in MES (their formula 38b and Figure 26a) for the same e. The correlation 
length ro depends in both calculation on the limiting magnitude. The redshift dependence 
of the CF is beside of the factor (1 + z)~ 3 ~ e also in ro- 

It is possible to concentrate the redshift dependence of the correlation function in the 
(1 + z) "i factor. By minimising the change of ro(m) with the limiting magnitude m, 
we get 

r 00 = 5.82 ± 0.05 and e = 2.60. (30) 



This result shows that the correlation length ro(z) decreases rapidly with the average 
redshift of the sub-catalogue. It means that the stable clustering is not yet reached, on 
small scales there is an inflow in physical coordinates. 

We estimate the average inflow/outflow velocities from the continuity equation as fol- 
lows: 

l-3^(^) 7 



1 + 



If r 3> ro v corresponds to the Hubble- flow and if r <C ro, then there is an inflow with 
v ~ —^^Hr. This makes possible e.g. to estimate the rate of galaxy merging. 
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4.5 The kernel of the integral @ 



In the following sections we determine at which r breaks the 3d CF down from the 
power law, by using the 2d CFs, measured from different sub-catalogues. To relate the 2d 
and 3d correlation functions of general form we must now determine the (||) kernel. The 
integral is easily calculated. The Q(t, m) = t 1_7 i^(t; m) functions are shown in Fig. 8 for 
15.5 < m < 18.5 mag. We included in Q(t, m) the multiplier t 1-7 , so the remaining factor 
under the integral is constant up to the point, where the 3d CF breaks down from scaling 
. With growing magnitude limit m, the maxima of the functions Q(t, m) shifted towards 
higher t = r9, they become flatter and in the integral (JsJ) the large t part of the 3d CF 
will be more important. 

The Q(t,m) functions have a very simple scaling. As it is to see from the non- 
relativistic, curvature free approximation of (||) 



Q(t,m) = bQ{t/a;m'). 

The coefficients for scaling on Q(t; 17mag) are given in Table 4. 

Table 4. Scaling factors to 17 mag 



(32) 



m 


a 


b 


loga 


log a 7 b 


15.5 


1.815 


0.187 


0.2589 


-0.2883 


16.0 


1.530 


0.302 


0.1847 


-0.2079 


16.5 


1.245 


0.545 


0.0952 


-0.1028 


17.0 


1.00 


1.00 


0. 


0. 


17.5 


0.844 


1.655 


-0.0731 


-0.0952 


18.0 


0.690 


3.141 


-0.1612 


0.2245 


18.5 


0.578 


5.467 


-0.2381 


0.3354 



Notice that b 



,-3 



The scaled functions are shown in Figure 9. By taking a 3d CF, 



independent of the sample depth, one can deduce from (32) a scaling law for the 2d CFs: 



w(9,m) = — w(a9;17). (33) 

In Figure 10. we show the 2d CFs measured from sub-catalogues m = 15.5 mag and 
m = 18.5 mag scaled to w(9; 17 mag) by using (32). One sees again that if the 3d CF 



independent from the depth of the catalogue, the large 9 part of the 2d CFs do not scale 
with the limiting magnitude. 

4.6 Modelling of the 3d correlation function 



Baugh 1996 solved the Limber-equation for £(r) by inverting (^|) through an iterative 
method. He used the full 17 < bj < 20 slice of the APM catalog. The resulting 3d 
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CF (for e 



4/;. 



-i 



< r < 



-3) well fitted by the (4.5Mpc/r) 1-7 power law for r < Ah 1 Mpc, between 
_ 25 Mpc £(r) has a shoulder and rising over the quoted power law. For 
r > 25/j- 1 Mpc the calculated 3d CF becomes negative and on the scales r > AOh 1 Mpc 
consistent with zero. Baugh 1996 calculates also the 3d CF for other magnitude slices and 
finds small discrepancies between these functions recovered from various magnitude slice 
of the APM catalogue. 

Here we use another method. We parametrise the 3d CF in different forms and compute 
the 2d CFs for several magnitude limits m from the Limber-equation (|8|) . The best fitting 
form and parameters determine the 3d correlation function. It is clear that the form of the 
3d CF must be the same type as that of the 2d CFs. To describe the physics mirrored by 
the correlation function, several distance like parameter are used: ro gives its amplitude 
for small r. This parameter is simple connected with the amplitude of the 2d correlation 
function (|24|). It is also clear that the 3d CF must cross zero at some r 2 , to permit the 
change of sign of the 2d CFs. It is generally assumed that there in no correlation between 
the galaxies at distances large enough ie for r — > oo the correlation function £ — > 0. This 
asymptotic is characterised by a parameter r^. Further, by considering some concrete 
form of the 3d CF £(r;ri,r3), one can not reconstruct with (|8|) the elbow of the 2d CF. 
We must introduce a further parameter in £ what describe the "shoulder" seen by Baugh 
96. We think, this shoulder corresponds to the extra correlation of the galaxies because 
of the cluster - cluster correlation. 

If we use the integral constraint / r 2 £(r)dr = 0, we get a relation between the param- 
eters. We have tried numerous function form. Here we describe only one parametrisation, 
since our main result, described in the next subsection do not depend sharply on the 
details of the parametrisation. Our best function has the form 



£(r) 



ro roY 1 _ fro n) v 
r T\) \r\ r 2 



-r/X 



(34) 



As it was shown in the previous section, one can not fit all the 2d CFs with a magnitude 
limit independent 3d CF. 

4.7 Three dimensional correlation function determined from different 2d sub-catalogues 



The parameters ro and e of the 3d CF are fixed by the scaling part of the 2d CFs and 
are given in (30). Now we determine the other parameters of the 3d CF £(r) in (34). For 
m= 17 mag we get n= 8.47 h _1 Mpc, r 2 = 39.95 h" 1 Mpc and A= 41.70 Ir 1 Mpc. The 
3d CF determined by these parameters is not far from the correlation function computed 
by Baugh 1996 from tha APM catalogue. However, if we use the same parameters for sub- 
catalogues with other magnitude limit, we find that for m < 17 mag the 2d CFs, calculated 
from (||) breaks down from scaling later than measured, for m > 17 mag earlier. We find 
that the 3d CF, calculated from different sub-catalogues, changes in a systematic manner. 
To stress that, we parametrise ri,r% and A as follows 



n = 7.50 + 0.65(m - 15.5)/i _i Mpc 
r 2 = 25.1 + 9.92(m - 15.5)/i~ 1 Mpc 
A = 31.2 + 7.98(m - 15.5)/i -1 Mpc 
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A = 13.3 + 5.10(m - 15.5)/i _1 Mpc 



(35) 



The 3d CF determined from subcatalogues with different limiting magnitude m are shown 
in Figure 11. The 3d CF determined by the parameters ri,rz and A do not satisfy the 
integral condition. To fulfil it, we should use in the last factor of (|34|) Ao instead of A, but 
it would give a bad fit for the negative part of the 2d CFs (the tail of the calculated 2d 
CF would be much less negative). It could be explained by the oscillation of the large r 
part of the 3d CF as in Tucker et al. 1997. 

The quality of the parametrisation and the choice of the parameters can be judged 
from the Figures 12 - 14. We have tried other dN/dz functions, other parametrisation for 
£, other e value (with depth dependent tq) and we found always that the 3d CF, calculated 
from the measured 2d CFs depend on the catalogue depth. 

5. Summary and discussion 

We determined the two point correlation function from the Muenster Red Sky Survey 
(MRSS, Ungruhe 1999). This survey uses the F-band and larger than the APM survey 
(MES). Maddox & al.1990 and MES have compared their measured correlation functions 
with those, determined from other catalogues. Here we compare our 2d CFs with the 
results of MES. 

The large dynamics range of the survey made it possible to determine the correlation 
function from subcatalogues with different magnitude limit m. The measured slope of 
the 2d CFs is somewhat larger for lower magnitude limit, but its average value 5 = 
7 — 1 = 0.69 ± 0.02, argees with the earlier measurements. The 9 max (m) value, where 
the CFs measured from subcatalogues with different magnitude limits m correspond to 
about 30/i -1 kpc projected distance, independent of m. When the projected distance of 
two galaxies is smaller than that, the software can not separate the two objects. The 
amplitude of the correlation functions, after the color correction, agrees excellently with 
those measured from the APM catalogue. As it is shown in Figure 3. our CFs agree well 
with those, measured by MES also in the region where the CFs do not follow the power 
law. Note that the break away from the power law begins at about the same angle, almost 
independent of the magnitude limit of the subcatalogue. After the zero crossing the MES 
correlation function is zero. We find that the amplitude of the negative part of the CFs is 
several times 10 -3 , and decreasing with increasing magnitude limit m (Fig. 12 - 14). The 
non- vanishing anticorrelation is a natural consequence of the definition of the correlation 
function. 

We have calculated from the measured 2d CFs the three dimensional two-point cor- 
relation function by using the Limber-equation. To determine the dN/dz function we 
used informations from 3d catalogues. We parametrizer it in the usual form and the 
parameters were for each subcatalogue deduced from the luminosity function of the Las 
Campanas survey (Lin et al. 1996). As it is shown in Fig. 9, the kernel of the integral in 

the w{9) £(r) transformation (deprojection) scales with the limiting magnitude of the 

subcatalogue. 

From a sub-catalogue with limiting magnitude m one can only determine a combination 
of ro(m) and e. There are several possibility to separate these parameters. The usual way 
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is to fix e on the basis of physical considerations or simple ignore it. To compare our 
parameters with the ro(m) and e, determined from the APM catalogue, we used e = 
— 1.3 and calculated r (m) from the subcatalogues. The agreement (Fig. 7) is excellent. 
However, if one considers e as evolution parameter it is more consequent to choose the 
correlation length tqq independent from the depth of the subcatalogue, ie. r (z) = r 00 (l + 

z) t With this parametrisation we have got roo = 5.82 ± 0.05/i Mpcande = 2.60. This 
result agrees with the general trend: deeper the catalogue lower the correlation length 
ro(z). This e value is unusually high. That means that there is yet an evolution in 
clustering. It is not easy to compare our ro,e with that of other authors, because the 
separation of the two parameters is not standardised. From the Las Campanas survey 
Jing, Mo & Borner 1998 found 7=1.86 and r =5.06 h" 1 Mpc. With e = 7 - 3 and z =0.1 
that gives ro ~ 4.6 Mpc, in agreement with our 5.82 • i.i - 5 - 6 / 1 - 68 — 4.2 Mpc. The limiting 
magnitude of the EPS catalogue (Guzzo & al., 1999) is bj < 19.4 mag; they found 7 = 
1.67 and r = 4.15 ± 0.2/i -1 Mpc. We get for r F < 18 mag r = 3.99 _1 Mpc. For deep 
surveys it is questionable if the parametrisation (||) correct. Interestingly, extrapolating 
our parametrisation ( |30[ ) to high z values we are not far from the measured values. E.g. 
z = 0.5 we get ro = 1.4 h _1 Mpc (in physical coordinates, go=0.5). This value agrees 
well with the result of Le Fevre & al.1996, (ro =1.33h~ 1 Mpc). Our roo,e values are 
compatible with the result of Postman & al. 1998, who find that if ro is fixed to be 
5.5±1.5 h -1 Mpc, then -0.4jej+1.3. Connolly & al. 1998 measure from a subcatalogue of 
HDF for fixed r = hAh~ l M V c e = 2.37^;^. 

By comparing our 3d CF determined from the r^jl8.5 mag catalogue with that of 
Baugh 1996 for 17 mag jfrjj20 mag, we see three minor differences. First, our correlation 
length ro (after taking into account the corrections for the different e's) is about a factor 
of 1.3 larger. It can be explained by the difference of the zq parameters of the dN/dz 
function. Second, our correlation function crosses zero at 50 h" 1 Mpc, that of Baugh at 40 
h _1 Mpc. Third, the CF of Baugh 1996 disappears at about 80 h _1 Mpc, ours seems to be 
negative also for larger distances. This is a consequence of the difference of the measured 
2d CFs in this region. 

As it is shown in Fig. 11 and ( |35| ) the position of the break-away from the scaling of 
the 3d CF and the zero-crossing varies with the depth of the sub-catalogues. This agrees 
with the result of Cappi et al. 1998. We conclude that the MRSS is not deep enough 
{z-med ~ 0.14) to determine up to which physical distance follows the 3d CF the power 
law. The structures, found in large and deep (Sloan and 2dF) catalogues are much larger 
then those seen in the MRSS. 

The author thanks all members of the Astronomical Institute of the University Muen- 
ster, especially Weltraud Seitter and Renko Ungruhe providing an excellent galaxy survey. 
I thank also George Efstathiou and Alex Szalay for valuable discussion. 
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Figure 1: Measured CFs for seven limiting magnitude from tf <15.5 to tf < 18.5 derived from 
152 plates 
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Figure 3: Correlation functions as measured from MRSS (small symboles) and from APM (lar; 
symboles) in MES Fig 23. 
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Figure 4: Comparison the APM and MRSS correlation function on linear scale 
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Figure 7: Correlation length ro in function of the limiting magnitude calculated from the APM 
catalogue and from MRSS for e—-1.3 
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Figure 8: The kernel of the integral (||) for different magnitude limits. 
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Figure 9: The function t 1 7 if (i; 17) and the other six t 1 7 K(t; m ^ 17) functions, scaled on it 
with @. 
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Figure 10: The 2CF, measured from a subsample < 17 mag and the 2CFs log w(9; 15.5) & 
log w(9] 18.5) scaled on it with (132" 
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Figure 11: 3d correlation function determined from the subcatalogues with limiting magnitudes 
15.5, 16.5, 17.5 and 18.5 mag. 
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Figure 13: The measured and and from the Limber-equation, using selection function 2 and 
in fl34| ) parametrised 3d CF calculated, 2d CFs for m—17.5 mag (upper panel), m=17 mag (lower 
panel) 
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Figure 14: The measured and and from the Limber-equation, using selection function 2 and 
in fl34| ) parametrised 3d CF calculated, 2d CFs for m=18.5 mag (upper panel), m=18 mag (lower 
panel) 
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